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Abstract: 


In  many  engineering  design  problems  it  is  possible  to  collect  data 
of  the  environmental  disturbances  which  are  acting  upon  our  systems. 

This  data  can  be  analyzed  by  determining  its  autocorrelation  and 
probability  density  function.  When  seeking  solutions  through  simulation 
it  is  desirable  to  be  able  to  generate  random  time  series  having  a 
predetermined  autocorrelation  and  probability  density.  This  paper 
describes  a  method  to  control  both  simultaneously. 

The  proposed  system  is  composed  of  a  linear  filter*  H{s) ,  with 
white  Gaussian  noise  as  input*  followed  by  a  nonlinear  element  £(x)* 
where  x(t)  is  the  output  of  H(s).  The  output  of  the  system  y(t)  is 
required  to  have  a  predetermined  probability  density  py(y)  and  a 
predetermined  normalized  autocorrelation  function  py(r).  The  non¬ 
linearity  f(x)  is  designed  to  give  the  required  density  py(y),  and 
is  relatively  easy  to  design  using  the  relationship 


f(x)  «  F  _1[F  (x) J 

y  * 


where  Fy  is  the  inverse  of  Fy(y)  which  is  the  cumulative  distribution 
function  of  y,  and 


FX(X) 


x  Jil 

f 


is  the  cumulative  Gaussian  distribution. 


j—  ...  , a  ..  _  _  .  ...  .  rn ,  .  .  .  . 


The  nonlinearity  f(x)  may  be  designed  manually  or  in  a  digital  computer 
by  using  a  double  table  look  up  or  in  a  hybrid  computer  system.  It  can 
then  be  stored  in  a  variable  diode  function  generator  of  an  analog 
computer. 

This  nonlinear  element,  however,  changes  the  autocorrelation  of 
the  input,  x(t).  The  amount  of  change  can  be  determined  by 


H  (x)  is  the  Hermite  polynomial  of  order  n. 
n 

This  calculation  can  be  done  most  conveniently  in  a  hybrid 
computer  system. 

The  normalized  autocorrelation  of  the  intermediate  stage,  x(t), 
px(r)  must  now  be  approximated  by  a  sum  of  exponentials: 

~slM  -s2M 

0x(r)  -  a^e  +  a2e  +  .... 

by  using  a  hybrid  computer  system.  One  approach  which  may  be  used  to 
determine  the  a^’s  and  s^'s  is  discussed  by  McDonough  amt  Huggins. 
The  filter  H(s)  may  be  obtained  from  the  a^*s  and  s^’s  by  spectral 
factorization. 


II 


By  forcing  the  system  with  white  Gai  ssian  noise  the  system 
output,  y(t)  will  have  the  desired  .autocorrelation  function  and 
probability  density. 
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1.  Introduction  j 


Most  of  our  physical  systems  are  subjected  to  environmental 
disturbances.  If  we  are  supposed  to  improve  our  system  through 
simulation  then  we  must  have  means  to  generate  disturbances  as  they 
appear  in  nature. 

In  most  cases  it  is  possible  to  record  these  random  disturbances 
for  further  analysis.  One  way  to  analyze  them  is  to  determine  their 
autocorrelation  function  and  their  probability  density.  Our  task  now  is 
to  design  a  system,  which  will  generate  a  random  time  series 
having  a  predetermined  autocorrelation  and  probability  density.  In  other 
words  we  must  be  able  to  control  both  simultaneously. 

Consider  a  system  described  in  Figure  1.  This  system  should  be 
so  designed  that  a  white  Gaussian  noise  as  an  input  will  give  y(t) 
as  an  output  having  the  desired  normalized  autocorrelation,  py(x)  and 
the  probability  density,  py(y).  Working  backwards  from  y(t),  the  instan¬ 
taneous  nonlinearity,  f(x)  must  be  such  as  to  guarantee  the  desired  proba¬ 
bility  density,  py(y) with  an  input,  x(t).  The  function  x(t)  is  always 
Gaussian,  but  it  must  have  a  normalized  autocorrelation,  Px(t)  such  as  to 
give  the  required  autocorrelation,  py(x)  after  going  through  the  nonlinear 
element  £(x).  The  linear  filter,  H(s)  must  be  so  designed  as  to  give 
Gaussian  output,  x(t)  with  zero  mean  and  the  required  autocorrelation  Px(x)  • 
For  simplicity,  the  variance  of  x(t)  will  be  normalized  to  unity.  The  non¬ 
linear  element,  f(x)  can  be  designed  first  because  p  (y)  and  p  (x)  are  known. 

y  x 

Then  P  (t)  can  be  determined  because  p  (x)  and  f(x)  are  known,  p  (x) 
x  y  x 

can  now  be  approximated  by  a  sum  of  exponentials.  And  finally,  the 
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Figure  1.  k  system  to  generate  random  time-series  having  a  predetermined 
autocorrelation  and  probability  density. 


linear  filter  may  be  obtained  from  the  a^'s  and  s^’s  by  spectral 
factorization.  We  will  now  have  a  system  which  will  continuously 
generate  random  time  series  having  the  required  properties. 

2.  Mathematical  Derivations. 

The  nonlinear  element  can  be  designed  by  using  the  relationship [3, 5, 6, 7] 
f(x)  -  F  ~X[F  (x)} 

y  x 


where: 

Fy  *  is  the  inverse  of  F  (y) 

Fy(y)  is  the  integral  of  the  desired  distribution  density,  p^(y). 

Fx(x)  is  the  integral  of  the  Gaussian  distribution  density,  p  (x)  or 
*  x 

X  -  £2 

V*>  -  L  j===  >*5 

V  2  TT 

Both  cumulative  density  functions  must  be  normalized  so  that  the  final 
values  of  the  integrals  are  one. 

The  relationship  of  the  normalized  autocorrelations  Px(t)  and  Py(t) 
can  be  determined  from  the  approximate  equation 
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Hn(x)  is  the  He  rial  to  polynomial  of  order  n[8].  The  number  of  terms  N,  is 
determined  hy  he  desired  normalized  error,  e,  in  the  approximation  which 


is  given  by 


N  K 


[i '  E  i2-  ] 


n=l  o^nl 


Since  p  (t)  and  f(x)  are  known,  Px(t)  can  be  derived.  To  design  the  filter, 
H(s)  the  normalized  autocorrelation  function  of  the  output,  Px(t)  must  be 
approximated  to  the  desired  level  of  accuracy.  This  can  be  done  by  finite 
auras  of  exponentials  [19]  as  described  by  R.  N.  McDonough  and  W.  H.  Huggins 
C 1 J -  In  their  approach,  "•••  a  givei-  number  n  of  exponentials  is  considered, 
such  that  the  integrated  squared  error  is  minimized  over  both  the  n  coeffi¬ 
cients  of  the  linear  combination  and  the  n  exponents  used.”  [from  the 
abstract  of  ref.  1.]. 

If  a  hybrid  computer,  however,  is  available,  then  the  mechanization  of 
finding  the  exponents  and  coefficients  can  be  simplified  as  discussed  later 
in  this  paper.  Also  it  is  believed  by  these  authors,  that  for  most  engineering 


cases  only  three  exponentials  are  required. 

The  problem  reduces  now  to  finding  the  best  estimate  of  p  (t),  P  (t)  of 

X  X 


the  form: 


-SjT  — S-T  -s_x 

Px(t)  =  a^e  +  a2e  +  a3e 


» 


These  three  exponents  can  be  easily  generated  in  an  analog  compitter 
and  by  using  optimizing  hybrid  techniques  the  parameters  a^'s  and  can 

be  so  adjusted  as  to  give  the  minimum  integrated  squared  error.  Then  the 
power  spectrum  would  be: 


X  2alSl  ,  2a2S2  ,  2a3S3 

$(s)  ®  2  2  +  22"^"  22 

-s  s2  -s  s^  -s 
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The  filter  H(s)  may  now  be  derived  by  spectral  factorization  which  for  the 
three  exponents  case  is  given  as  follows:  [2,  4] 


As4  -  2Bs2  +  C 


„2  „  B  ±  VB^  -  AC 
A 


B  +  V  o  -  AC 


U2  -  AC 


where 


A  «  2(a^s^  +  &2s2  +  a3s3^ 

2  2  2  2  ?  2 
B  **  a^s^s^  +  +  a2S2s^  +  826282  +  a3s3si  +  323382 

2  2  2  2  9  2 

C  *  2 (a^sis2  ®3  +  a2s2sl  s3  +  a3®3sl  ®2  ) 

Wfswt\/2alsl  +  2a2s2  +  2a3s3l<s  +  Pl>  (s  +  P2> 

(S  +  Sx)  (S  +  S2)  (S  +  S3) 

or  more  conveniently,  defining 


y2a^s^  +  2&2s2  ***  ^a3s3 

Q  (Pi~sl>  (pg-si)  +  a(p1-s2)  (P2~s2)  o(P1-s3)(p2-s3) 

(s+sj) (s2”S^) (s3-Sj)  (s+s2) (s^— s2) <s3-s2)  (s+s3) ^si“s3) (s2”s3) 


In  the  general  case  with  M  exponents,  we  have 


p  <t)  »  Z  a.  e"SiT 
x  i-1  1 


Then  the  power  spectrum  may  be  written  in  the  form 
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v*  '-*---4  'v  r  ¥xr  ‘A 


M-l  „ 

M  2a  s  /  (pi  "  8  * 

4>(s)  =  E  _ 1~-1--  g  M _ 

vls;  .  .  ,  2  2.  M  „  0 

i=l  '8.  —  s  )  ,2  2. 

i  n  (S.  -  s  ) 

i-1  1 


and  therefore,  the  filter  is  given  by 


H(s) 


M-l 

IT  (p.,  +  s) 
i=l 
M 

ji  (s.  +  s) 
i=l 


It  is  easy  to  see,  that  to  build  the  filter  the  same  analog  computer 
circuit,  which  was  used  to  find  the  a^’s  and  s^'s,  can  be  used  just  by 

changing  the  setting  of  the  three  input  potentiometers. 

To  clarify  the  method  an  example  with  an  analytical  solution  is  given 
in  Appendix  No.  1. 


3.  Computerized  Methods. 


When  seeking  solutions  for  these  types  of  problems  through  hybrid 
computation  the  following  steps  have  to  be  taken.  First  of  all,  two 
distinctively  different  stages  can  be  defined,  namely:  the  design 
stage  and  the  stage  of  utilization  of  the  designed  system  to  generate 
random  time  series. 

A.  Design  Stage. 

1)  The  desired  probability  density  function  p^(y)  is  given 
either  in  a  numerical  or  in  a  graph  form.  In  the  first 
case  quite  often  a  digital  computer  subroutine  is 
available.  In  many  cases  it  can  be  conveniently 
generated  in  an  analog  computer.  In  the  second  case  a 
curve  follower  or  magnetic  Ink  could  be  used.  In  any 
event,  the  given  Py(y)  must  be  generated,  digitalized, 
stored  and  plotted.  It  must  be  repetitively  available 
for  further  study.  6 


2)  Gaussian  probability  density  function  in  scat  cases  is 
available  as  a  digital  computer  subroutine.  It  must  be 
generated  and  stored  as  in  No.  1. 

3)  The  stored  p^(y)  must  be  fed  into  analog  computer  for 
integration  to  form  F^(y).  This  will  then  be  sampled  and 
stored  in  a  digital  computer. 

4)  Gaussian  density  function,  px(x)  is  integrated  and  stored 

as  in  No.  3  to  form  F  (x). 

x 

5)  The  given  autocorrelation  function  of  the  output,  y(t)  in 

its  normalized  form  p^(t)  must  be  generated,  digitalized, 

stored  and  plotted.  Notice  step  No.  1. 

in)  Because  F  (y)  and  F  (x)  now  are  available  the  non-linear  ele- 
y  x 

ment,  f(x)  can  be  designed  and  stored  in  a  digital  computer. 
This  can  be  done  by  using  a  double-look-up  method  in  a  digital 
computer  as  Illustrated  in  Figure  2.  In  the  final  run*  either 
a  digital  computer  or  a  variable  diode  function  generator  in  an 
analog  computer  could  be  used  to  implement  the  required  non¬ 
linear  element,  f(x). 

7)  To  generate  K^'s,  n  •  0,1, ...,N,  we  need  f(x),  p^(x)  and  the 
Hermite  polynomial,  (x)  .  This  can  be  done  most  conveniently 
in  an  analog  computer  as  of  Figure  3.  Notice  that  two  previous 
stages  of  fi(x)  must  be  stored  for  the  next  set  of  K's.  This 
should  be  continued  until  the  error  e  is  satisfactorily  small. 
[Hermite  polynomials  -  H^(x)  *  xHq (x)  -  nHn_^(x);  Kq  *  1, 

Hj  *  x,  Hj  *  -  1,  -  3x.]  Kq  is  left  in  this  circuit 

for  checking  purposes,  since  Kq  is  equal  to  the  mean  of  y,  i.e. 

OB 

Kom  f  yVy)dy- 

—w  * 
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Figure  3. 


An  analog  coaputer  circuit  to  generate  K. 

x2 

*n  -  7  £<x)  ^(x)  dx 

—  ^2* 

£(x)  stored  in  digital  line  1. 


-  x_ 
e  2 


stored  in  digital  line  2. 


8)  The  normalized  autocorrelation  function  of  x(t),  Px(t) 
can  now  be  generated,  digitalized,  stored  and  plotted. 

An  illustration  of  how  this  can  be  done  in  a  case 

where  four  values  of  were  found  to  •satisfactory,  see 
Figure  4.  pxCr)  can  be  generated  in  C  by  letting  a 
servo-mult iplyer.  A,  maintain  the  output,  B,  at  zero. 

9)  To  approximate  px(t)  by  a  sum  of  exponentials  the 
following  method  is  used  as  a  first  trial.  Build 
exponential  generating  circuits  in  an  analog  computer, 
for  example  three  of  them.  Take  the  absolute  value  of 
error  between  the  sum  of  the  exponents  and  the  stored  value 
of  Px(t)  and  integrate  in  an  analog  computer.  By  using 
digital  computer  optimizing  techniques  readjust  the 
servoset  potentiometers  representing  the  a^s  and  s^'s  as 
to  minimize  the  error.  Print  out  the  optimizing  values 

of  a^'s  and  s^'s. 

10)  Design  the  filter  H(s)  as  per  Chapter  2.  Notice  that  the 
analog  computer  circuit  in  Step  9  can  be  used  by  merely 
changing  the  values  of  the  input  potentiometers  (and  of 
course  by  removing  the  initial  conditions) . 

B.  Stage  of  Generation  of  Random  Time  Series 

1)  Apply  white  Gaussian  noise  to  the  system.  Plot  f>y(t)  and 
compare  with  the  given  Py(t). 

2)  Plot  Py(y)  and  compare  with  the  given  Py(y). 

3)  Plot  the  impulse  response  of  the  filter,  H(s). 

4)  Plot  px(r). 
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4.  Experimental  Computer  Run. 

For  this  experimentation  an  example,  which  has  been  analytically 
solved  in  Appendix  1,,  was  chosen.  This  example,  however,  was  originally 
designed  for  a  fast  hybrid  system  with  a  large  memory.  To  fit  the  system 
available,  the  hvbrid  computer  system  of  «-he  Department  of  Industrial  and 
Systems  Engineering;  University  of  Florida,  Gainesville,  Florida  [10.] 
the  time  constant  had  to  be  increased.  The  required  autocorrelation  and 
probability  density  of  the  output,  y(t)  were 

P^(t)  »  e  where  6  ".5 

P  (y)  ®  uniform  between  -  1  and  +  1. 

The  linear  filter,  H(s)  changes  now  to 


H(s) 


1.06 
s+.  5 


.119 

s+1.5 


and  Figure  5. a  shows,  how  the  non-linear  element,  f(x)  was  designed*  This 
was  stored  in  the  variable  diode  function  generator.  A  permanently  stored 
subroutine  war  used  to  generate  random-numbers  having  a  uniform  distribution 
[Fig.  5.b.).  Sequentially  eight  of  these  were  averaged  to  give  one  value 
for  the  required  normally  distributed  data  [Fig.  5.c» ].  As  one  can  easily 
see  there  is  a  long  range  fluctuation  present  in  this  data.  This  irregular 
fluctuation  could  be  called  a  DC-component  and  its  effect  is  detrimental 
for  the  correllogram  when  one  is  forced  to  use  such  a  short  record 
of  data,  only  630  data  points,  which  means  that  the  base  of  cn  autocor¬ 
relation  function  can  be  only  315  data  points.  To  remove  at  least  part 
of  this  DC-compoment  a  double,  exponential  smoothing  was  used.  This  means 
that  the  data  [Fig.  5.c.}  was  run  through  an  exponential  smoothing  circuit 
with  a  time  constant  of  1  sec.,  then  the  smoothed  data  was  run  again 
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Figure  5.b 


backwards  through  the  same  circuit  to  compensate  for  the  shift.  This 
doubly  smoothed  data  was  then  subtracted  from  the  original  Gaussian  data 
to  form  the  final  white  noise,  w(t)  to  be  used  as  the  input  to  the  linear 
filter,  H(F)  (Fig.  5.d.).  It  is  better,  though  not  perfect.  The  output 
of  the  filter,  x(t)  (Fig.  5.c.]  was  then  adjusted  in  such  a  way  as  to  have 
a  standard  deviation  of  25  volts  and  its  autocorrelation  function  was 
plotted  [Fig.  5 . f . } .  x(t)  was  then  fed  into  the  non-linear  element,  f(x) 
to  give  an  output,  v(t)  [Fig.  5.g.].  Hie.  autocorrelation  function,  Py("0 
was  determined  [Fig.  5.h,]  and  the  histogram  of  y(t)  plotted  [Fig.  5.i.]. 
The  hybrid  computer  system  used  in  this  experimentation  and  the  method  for 
obtaining  the  autocorrollographs  have  been  described  earlier  [10],  [11]. 

Considering  the  limited  number  of  data  points,  which  could  be  used 
in  tliis  hybrid  system,  the  histogram  is  surprisingly  good.  The  autocorre¬ 
lation  plots  show  that  the  removal  of  the  DC-component  was  not  as  good  as 
was  hoped  for.  However,  this  discrepancy  would  disappear  in  time  if  a 
continuous  Gaussian  noise  source  were  available  and  longer  record  could  be 
used  to  check  the  autocorrelation  functions. 

5.  Concluding  Remarks. 

A  method  has  been  developed  to  generate  random  time  series  having  a 
predetermined  autocorrelation  and  probability  density.  Needless  to  say 
that  this  research  is  still  in  its  infancy.  Also  due  to  the  lack  of  space 
wc  are  unable  to  show  some  interesting  results  concerning  the  shape  of  the 
non-linearities  and  their  effects  to  the  autocorrelation  function  of  x(t). 
This  research  will  definitely  continue  especially  to  find  out  what  are  the 
limitations,  if  any,  of  this  approach. 
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Figure  5.d 


ttput  of  the  linear  filter  B(s),  x(t),  having  a  noma: 
predetermined  autocorrelation  function  p  (c). 


The  histogram  of  y(t)  showing  the  computed  distribution  py(y)  as  compared  to  the 
required  uniform  distribution. 

’  '  ’*  .  I 


Appendix  No.  1 


A  ’Simple  Example. 


An  example,  where  an  analytical  solution  can  be  found  is  given 
to  test  the  validity  of  this  method.  Suppose  that  the  output,  y(t)  is 
required  to  have  a  probability  density  function  of  the  form:  [8] 

P  (y)  »  §£  exp  (a2^l)[erf“‘i(|^)]2}  ,  I  yj£A 

'  *  \ 

where  erf  *  is  the  inverse  cf  the  function  ;■  , 


■  x  ~  2  '-A';-  A-  .  - 

er€<x)  «  /  --T— >-  du 

;;  ::  o  -  V*.  '  .  _  r-  -•/“* 

Let  the  desired  autocorrelation  function  of  y(tj  be  of  the  form 

,  i  =  c  ~B.lt  I  y  ,  / 

;  Ry(x)  «  Ry(0)  e  , 


or  in  normalized  form: 

-b|t| 

Py(x)  -  e 


In  this  case  the  variance  of  y(t)  is 


o2  -  R  (0)  »  ~  arc  sin  — «=• 
y  J  *  l*az 


For  a  «  1  the  density  of  y(t)  is  uniform  in  (~A,A) .  - 


2A2  Px<T)  . 

*y(v  ■  —  “rc  8l” 


White,  noir.e  with  R^,(t)  =  cS(t)  is  used  a?  the  input  to  the  tilter  and 
the  non-linear  element  *■(•')  is  of  the  form. 

f(>;)  -  2A  fXt  (-~) 

where  u2 

erf  (x)  =  /  e  .  2  -  ■  3« 

0  VfiT  .  \ 

Numerical  example  follows.  -  _  '  ..  .  ;  ; 

Let  us  assume  that:  -  .  : 

a  «  1  and  8  «=  5  then  ,  "  -  :  ~  .  .  '  '  .  .  '  .  - 

K  «  arc  sin  — *•  7-'*?  *524  ‘ :V  >  . '  -  ,  --V  ' 

1+ci2  6  -y:  ^  V \ ; :.v  . ',:v  • 

K,  =  73.6;  K3  *  4.56  and  K2  =  37*5  •  Y Y>\  ^ 

B  *  5.0C  B  **  *575  and  B  «  .0125  •  ■  ,  -  ' ;  .  .. 

1  2  J  -r  :  .  -  :  . 

f5.08  .575  .0125 '  . 

HCt)  “  .66^s+5  -  s+i5  +  s+25j 

3.34  ,  .378  +  .0082 

s+5  s+15  s+25 


Only  two  terras  required  in  this  case 


Co.'ifjiX,  io ni.  1  y 


v  (i)  "•  (I  +  «“)  «J.n  • — <’- 


wK ,(■'  ) 
2 


2  A 


vh  ich  is  a pern:: Lina  C ely : 


P  (t> 
x 


(i  +  &2)\ - „  K  (t) 

(2A“  ^ 

L 


r  it 


?.A 


2*\ 

y 


R..CO 


-to  r  b 

‘T  .  i; 


2  VT>f  ) 

2A  *  J 


3! 


5! 


These  terms  provide  enough  accuracy. 
Then  if  we  let: 


K  (0)n 

y 


=  K  =  arc  sin 


2P/ 


1+ct 


0  <  K  < 


N  *“ J 


we  have: 


PX(T> 


,  r-eWi  k2  -33M  k«  -Mhl? 

(1  +  a“)  K  |e  -  TT-  e  Hgg-  *  J 


K 


K 


i.xcs)  -  a  +  «o  Kf.  +  ggp 


v' 


_  (j  +  a2)  Kg  fe(9g2-s2)(25g2-92)  -  12k2(32-s2)  (2582-b2) 
12  L  (g2-s2)  (9g2-s2)  (25g2-s2) 

+  (S2-s2)  (962-s2)I 
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,0.  -to2)  Kg  :g^|24-mWj  -  2sV-  im  -  1S6K?-  4 
.12(62- s2)  (9jJ2-s2) 

5  A  •'  15400  "  30QK2  i-  9K4T  :V  j 


This  vili  result  in  the  required  filter  as: 


=  fcT-f  £x  m  . 


V  12 


(3+s)  (30  is)  (50+s) 


where 


K,  =  >y/  5400  -  300  1C2  +  9  K4,  .  .real  tor  0<K£  ~ 


K„  «  \f  24  -  12  K2  +  K4,  real  for  .0  £  K  <  ^ 


K2  *  v/2  { (40S-15G  -K?  +  3  K4)  +  Kx  K^},  real  for  0  <  K 


U(s)  =  yj 


^(i+a2)  K2 

r  o 

l 

12  '•  i 

1  f+3 
^  - 

where 


,  ..  V  S  +  S 
V~~ 8 - 


Kl“  3K2  +  9K* 
- 


K  -  5k2  t  25K3 


•f  If v  *■? i *:  *  ’ v y gvwto}  V'*’ 
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1J.  ABSTRACT 


^  In  many  engineering  design  problems  it  is  possible,  to  collect* data  of  the 

**  *  „  j'  j»  '  '  - ‘V_  ’  ' 

environmental  disturbances  which  are  acting  upon  our*  systems.  This  data  can  be  ! 
analyzed  by  determining  its  autocorrelation  and  probability  density  function. 

Uhen  seeking  solutions  through  simulation  it  is  desirable  to  be  able  to  generate 
random  time  series  having  a  predetermined  autocorrelation  and  probability  density. 

€L  - 

This  paper  describes  a  method  to  control  both  simultaneously.  /  )  isar- 

The  proposed  system  is  composed  of  a  linear  filter,  H(s),  with  white  Gaussian 
noise  as  input,  followed  by  a  nonlinear  element  £(2),  where  x(t)  is  the  output  of 
H(s).  The  output  of  the  system  y(t)  is  required  to  have  a  predetermined  probability 
density  py(y)  and  a  predetermined  normalized  autocorrelation  runction  p  (t).  The 
nonlinearity  £(x)  is  designed  to  give  the  required  density  py(y),  and  is  relatively 
easy  to  design  using  the  relationship 

«x5  -  Fy'V.frH 

DD  .iSK. 1473 
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Abstract  (continued) 

where  l7  is  the  Inverse  of  F  (y)  which  is  the  cumulative  distribution  function 

v  y 

of  y,  and 

x  - 

f;<X>  "  S7==^ 

~ “  v2* 


is  the  cumulative  Gaussian  distribution. 

The  nonlinearity  £(x)  may  be  designed  manually  or  in  a  digital  computer  by  using 
a  double  table  look  up  or  in  a  hybrid  computer  system.  It  can  then  be  stored  in 
a  variable  diode  function  generator  of  an  analog  computer. 

This  nonlinear  element,  however,  changes  the  autocorrelation  of  the  input, 
x(t) .  The  amount  of  change  can  be  determined  by 

K  K  ^ 

pv(t)  -  i  -r-  Ip  (t)]n 
y  n»l  n\  x 


H  (x)  is  the  flermite  polynomial  of  order  n.. 

This  calculation  can  be  done  most  conveniently  in  a  hybrid  computer  system. 
The  normalised  autocorrelation  of  the  intermediate  stage,  x(t),  p^Ct)  must 
now  be  approximated  by  a  rna  of  exponentials:  : 

*  -Sl{Tj  S2\*\ 

Pjj(^)  *  +  a2e  +  ... 

by  using  a  hybrid  computer  system.  One  approach  which  may  be  used  to  determine 
the  a^'e  and  s^’s  is  discussed  by  McDonough  and  Huggins.  The  filter  H(s)  may  be 
obtained  from  the  a^s  and  s^s  by  spectral  factorization. 

By  forcing  the  system  with  white  Gaussian  noise  the  system  output,  y(t)  will 
have  the  desired  autocorrelation  function  and  probability  density. 


